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Among the major uncertainties involved in the Chandrasekhar mass models for Type 
la supernovae are the companion star of the accreting white dwarf (or the accretion rate 
that determines the carbon ignition density) and the flame speed after ignition. We 
present nucleosynthesis results from relatively slow deflagration (1.5 - 3 % of the sound 
speed) to constrain the rate of accretion from the companion star. Because of electron 
capture, a signihcant amount of neutron-rich species such as ®‘^Cr, ^’’Ti, ®®Fe, ®^Ni, etc. are 
synthesized in the central region. To avoid the too large ratios of ^^Cr/^®Fe and ^’’Ti/^^Fe, 
the central density of the white dwarf at thermonuclear runaway must be as low as C 2 
X 10® g cm“^. Such a low central density can be realized by the accretion as fast as 
M ^ 1 X 10“’’ Mq yr“^. These rapidly accreting white dwarfs might correspond to the 
super-soft X-ray sources. 

1. INTRODUCTION 

Supernovae of different types have different progenitors, thus producing different heavy 
elements on different time scales during the chemical evolution of galaxies. Because the 
lifetime of their massive progenitors is about 10®“’’ years being much shorter than the age 
of galaxies. Type 11 supernovae (SNe 11) and Type Ib/lc supernovae (SNe Ib/lc) cause 
the heavy-element enrichment in the early phase of the galactic evolution. 

In contrast. Type la supernovae (SNe la) produce heavy elements on a much longer 
time scale in the later phase of the galactic evolution. There are strong observational 
and theoretical indications that SNe la are thermonuclear explosions of accreting white 
dwarfs (e.g., Nomoto et al. 1994 for a review). However, the exact binary evolution 
that leads to SNe la has not been identihed (e.g., Renzini 1996; Branch et al. 1995 for 
recent reviews). The identihcation of the progenitor’s evolution is critically important 
1) for clarifying whether the nature of SNe la at high redshift is the same as nearby 
SNe la, and 2) for understanding the origin of some systematic variations of light curves 
(i.e., brighter-slower), which is a critical material for the determination of cosmological 
parameters, Hq and Qq. In this paper, we provide some important constraints on the 
progenitor system from the viewpoint of nucleosynthesis, namely, the carbon ignition 
density which is translated into the accretion rate for the Chandrasekhar mass models. 
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■ skhar iffass whitp^warf model, carbon burning in the central region 
cleai/runwajyA flame front then propagates at a subsonic speed as a 
1^0 heat/tfansport across the front (e.g., Nomoto et al. 1996a for a 
fjor uncertainty is the flame speed which depends on the development 
xioHS scales at the flame front. Multi-dimensional hydrodynamical 
ime propagation have been attempted by several groups (Livne 1993; 
4; Khokh^Q^95; Niemeyer & Hillebrandt 1995). These simulations 
a carbon deflagration wave propagates at a speed Udef as slow as 1.5 - 
;ed Vs in the central region of the white dwarf. Though the calculated 
flame ^eed is still very preliminary, it is useful to examine nucleosynthesis consequences 
of si^ a slow flame speed. 

I^the deflagration wave, electron capture enhances neutron excess, which depends on 
Kom Udef and the central density of the white dwarf pg = Pc/10® g cm“^. The resultant 
nucleosynthesis in slow deflagration has some distinct features compared with the faster 
peflagration like W7 (Nomoto et al. 1984; Thielemann et al. 1986), thus providing 
impg)gf 0 ]:|^constraints on these two parameters. The constraint on the central density is 
equivalent to a constraint on the accretion rate. 

We calculate explosivei nucleosynthesis for two cases with po = 1.37 (C) and 2.121 (W) 
^at the thermonuclear runaway, i.e., at^t^e stage when the timescale of the temp^ature 
"^rise in the center becomes shorter th^' the dynamical tirnj^^^ Here C and W imply 
that the models are the same as calculated for C 6 and WT; respecgtively (Nomoto et al. 
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Figure 2. The abundance distribution in slow deflagration WS15. 


1984). For the slow (S) deflagration, we adopt v^ef/vs = 0.015 (WS15, CS15) and 0.03 
(WS30). The central region behind the slower deflagration undergoes electron capture for 
a longer time than in W7, thereby having significantly reduced Yg. In Figure |l|, profiles 
of Ye for these cases and W7 are shown (Thielemann et al. 1996b). In general it can be 
recognized that small burning front velocities lead to steep 1^-gradient which flatten with 
increasing velocities. Lower central ignition densities shift the curves up (CS15), but the 
gradient is the same for the same propagation speed. 

Figure ^ shows the abundance distribution of neutron-rich species such as ®^Cr, ^°Ti, 
^®Fe, and ®^Ni behind the slow deflagration for WS15. The locations of ®‘^Fe and ®®Ni, 
overproduced in W7, correspond to Yg = 0.47-0.485. Due to the De-gradients which are 
steeper than for W7, the amount of matter in a given De-range is reduced, but also smaller 
central values are attained, giving rise to more neutron-rich nuclei. Yg ~ 0.46 ~ 26/56 
leads to a large abundance of stable ®®Fe (not from ^®Ni decay). Yg = 0.44 - 0.46 result also 
in '^^Ca, ^°Ti, ^'^Cr, and ®®Fe. Of these nuclei '^^Ca with Z/A 0.42 is only produced if De 
approaches values close to and smaller than 0.44 (Woosley 1996; Meyer et al. 1996). As 
the deflagration wave propagates outwards, the white dwarf gradually expands to undergo 
less electron capture and thus mostly ^®Ni is synthesized. Eventually, the deflagration 
enters the region of incomplete Si burning and explosive 0-Ne-C burning. 

The neutron excess is sensitive to the initial central density. For WS15 (pg = 2.12), the 
masses of ^^Cr and ^°Ti are smaller than those for CS15 (pg = 1.37) by a factor of ~ 5 and 
20, respectively, and those for WS30 (vdet/vs = 0.03) by a factor of 2 and 4, respectively. 
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Figure 3. The abundance distribution in the delayed detonation model WDD2 as a 
function of interior mass and the expansion velocity. 


3. NUCLEOSYNTHESIS IN DELAYED DETONATION 

If the deflagration speed continues to be much slower than in W7, the white dwarf 
undergoes strong pulsation as hrst found by Nomoto et al. (1976). In this pulsating 
deflagration model, the white dwarf expands to quench nuclear burning when the total 
energy of the star is still negative. Then the star contracts to burn more material to make 
the total energy positive (~ 5 x 1049 gj-g g 

The deflagration might induce a detonation at low density layers. In the delayed detona¬ 
tion model (Khokhlov 1991a; Woosley & Weaver 1994), the deflagration wave is assumed 
to be transformed into detonation at a certain layer during the hrst expansion phase. In 
the pulsating delayed detonation model (Khokhlov 1991b), the transition into detonation 
is assumed to occur near the maximum compression due to mixing. 

We study explosive nucleosynthesis with a large reaction network (Thielemann et al. 
1996) for various delayed detonation (DD) models, which has not been discussed in detail 
before (Nomoto et al. 1996b). We artihcially transform the slow dehagration WS15 into 
detonation when the density ahead of the hame decreases to 3.0, 2.2, and 1.7 x 10^ g 
cm“^ (WDD3, WDD2, and WDDl, respectively, where 3, 2, and 1 indicate pj at the 
transition). Then the carbon detonation propagates through the layers with p < 10® g 
cm“®. The explosion energy of three WDD models is 1.5 x 10®^ ergs s“^ and the mass 
of '^^Ni is 0.73 Mq (WDD3), 0.58 Mq (WDD2), and 0.45 Mq (WDDl). 

Figures |] shows the abundance distribution against the expansion velocity and Mr after 
the passage of the slow deflagration (WS15) and the delayed detonation. For comparison, 
the abundance distribution of W70 is shown in Figure H, where the initial metallicity (i.e., 
the initial CNO elements which are later transformed into ^^Ne) is assumed to be zero. 
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Figure 4. The abundance distribution in the delayed detonation model WTO as a function 
of interior mass and the expansion velocity. 


It is seen that WDD2 and WDDl produce two Si-S-Ar peaks at low velocity (~ 4000 
km s“^) and high velocities (10,000 - 15,000 km s“^). The low velocity intermediate mass 
elements are important to observe at late times to distinguish models. In particular, the 
minimum velocity of Ca in WDD models is ~ 4000 km s“^, which should be compared 
with the observed minimum velocities of Ca indicated by the red edge of the Ca II H and 
K absorption blend (Fisher et al. 1995). 

Meikle et al. (1996) have observed a P Cyg-like feature at ~ 1.05/1.08 pm in SN 1994D 
and 1991T. They note that, if this feature is due to He, He in SN 1994D is likely to be 
formed in a-rich freezeout and mixed out to the high velocity layers (~ 12,000 km s“^). 
The maximum velocity of He is 5000 - 6000 km s“^ in WDDs being slower than ~ 9000 
km s“^ in W7, so that more extensive mixing of He would be required for WDDs than 
W7. Alternatively, if the feature is due to Mg, the Mg velocity is conhned in 12,500 - 
16,000 km s“^ in SN 1994D, which is consistent with W7 (13,000 - 15,000 km s“^). For 
WDDs, on the other hand, the minimum velocities of Mg are 14,500 km s“^ (WDDl), 
16,500 km s“^ (WDD2), and 18,000 km s“^ (WDDS), and the latter two models seem 
to have too high velocities. 

4. YIELDS OF TYPE lA SUPERNOVAE 

Total isotopic compositions of WDDl - WDDS are given in Table |^, and compared 
with the solar abundances in Figure ^ which are normalized to ®®Fe. Table |I| also includes 
W7 and WTO updated with the latest reaction rates along with WTO (Thielemann et al. 
1996b). We note: 

1) The synthesized amounts of Fe and thus the ratio between the intermediate mass 
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Figure 5. The ratios of integrated abur^ances 
after decay of unstable nuclei, normalize^ tOj ^ 
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elements to Fe, Si-Ca/Fe, are sensitive to ^lie transition oensi 
detonation. Among the WDD models, WDp2 produces almost 
^®Ni as W7 (~ 0.6 M©) but more Si-Ca than W7 by a factor of 
oxygen is burned in the outer layers. WDDl fiiis even larger Sh 
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to the solar ratio. Therefore, the measuremenir' 

(Tycho, SN 1006, etc.) would be useful to distin, 

2 ) Some neutron-rich species such as and ^"^Ti are overproduced with respect to 
®®Fe. To see the degree of overproduction, we combine nucleosynthesis products of SNe 
la and SNe II with various ratios and compare with solar abundances of heavy elements 
and their isotopes. 

Nucleosynthesis products of SNe II as a function of stellar masses are taken from the 
calculations by Nomoto & Hashimoto (1988), Hashimoto et al. (1989, 1996) and Thiele¬ 
mann et al. (1996a) as summarized in Tsujimoto et al. (1995). SNe II yields integrated 
over mi = 10 Mq to = 50 Mq with the Salpeter IMF are given in Table |^. The upper 
mass bound is chosen to give [0/Fe] = -1-0.4, which is consistent with the observations 
of low metallicity stars for [Fe/H] < -1. 

Nucleosynthesis products of SNe la are those from WDD2 and W70, and the results 
are shown in Figures l-ll- For WDD2, the best £t to the solar abundances are obtained 
for the ratio between SNe la and SNe II contributions as r = 0.07 where r is the mass 
fraction contributed by SNe la per unit mass of all heavy elements in the gas as dehned 
in Tsujimoto et al. (1995). 

Aided with a reasonable chemical evolution model (Yoshii et al. 1996), the number of 
SNe la ever occurred relative to SNe II is determined to be Nh^/Nu = 0.12 in order to 
reproduce the observed abundances. This is consistent with the fact that the observed 
estimate of SNe la frequency is as low as 10 % of the total supernova occurrence (van den 
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Figure 6. Solar abundance pattern based on synthesized heavy elements from a com¬ 
posite of Type la and Type II supernova explosions with the most probable ratio of r. 
Relative abundances of synthesized heavy elements and their isotopes normalized to the 
corresponding solar abundances are shown by circles. Here WDD2 is adopted for the 
Type la supernova model. 



Figure 7. Same as Figure | but for W70 as an SN la model. 



















Bergh & Tammann 1991; Tsujimoto et al. 1995). 

With this relative frequency, ^®Fe from SNe la is about 50 % of total ®®Fe. Then the 
abundance ratios between neutron-rich species and ^®Fe could be reduced by a factor of ~ 

2. For W70, the excess of ®®Ni/®®Fe is now within the uncertainties (Fig. |^). For WDDs, 
on the other hand, ^^Cr and ®°Ti are overproduced as seen in Figure ^ even when the 
contribution from SNe II is taken into account. 

The above results imply that the central density of the Chandrasekhar mass white 
dwarf at thermonuclear runaway must be as low as 2 x 10® g cm“^, though the exact 
constraint depends somewhat on the flame speed. Such a low central density can be 
realized by the accretion as fast as M > 1 x 10 “'^Mq yr“^. Such rapidly accreting white 
dwarfs might correspond to the super-soft X-ray sources. 
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Table 1 

Nucleosynthesis products of SN 11 and la 


Species 


Synthesized 

mass (Mq) 



Type II 
1O~5OM0 

W70 

W7 

Type la 
WDDl 

WDD2 

WDD3 

-TIq 

7.93E-02 

5.08E-02 

4.83E-02 

1.80E-03 

1.15E-03 

5.82E-04 

13c 

3.80E-09 

1.56E-09 

1.40E-06 

3.52E-08 

1.04E-08 

6.90E-09 

14n 

1.56E-03 

3.31E-08 

1.16E-06 

3.08E-04 

2.44E-04 

7.42E-05 

15n 

1.66E-08 

4.13E-07 

1.32E-09 

6.84E-07 

2.75E-07 

1.68E-07 

160 

1.80 

1.33E-01 

1.43E-01 

9.96E-02 

6.93E-02 

4.69E-02 

1^0 

9.88E-08 

3.33E-10 

3.54E-08 

3.99E-06 

4.50E-06 

1.72E-06 

180 

4.6IE-03 

2.69E-10 

8.25E-10 

6.96E-07 

4.62E-07 

2.42E-07 

19p 

1.16E-09 

1.37E-10 

5.67E-10 

1.68E-09 

8.90E-10 

5.48E-10 

20Ne 

2.12E-01 

2.29E-03 

2.02E-03 

1.45E-03 

9.13E-04 

6.77E-04 

2iNe 

1.08E-03 

2.81E-08 

8.46E-06 

4.09E-06 

1.47E-06 

2.30E-06 

22Ne 

1.83E-02 

2.15E-08 

2.49E-03 

1.34E-05 

1.96E-06 

1.39E-06 

to 

CO 

6.51E-03 

1.41E-05 

6.32E-05 

3.20E-05 

1.30E-05 

6.53E-06 

24 Mg 

8.83E-02 

1.58E-02 

8.50E-03 

8.29E-03 

4.76E-03 

2.93E-03 

25 Mg 

1.44E-02 

1.64E-07 

4.05E-05 

4.60E-05 

2.39E-05 

1.44E-05 

26 Mg 

2.0IE-02 

1.87E-07 

3.18E-05 

5.52E-05 

3.57E-05 

1.09E-05 

27 Al 

1.48E-02 

1.13E-04 

9.86E-04 

4.65E-04 

2.74E-04 

1.73E-04 

28 Si 

1.05E-01 

1.38E-01 

1.50E-01 

3.48E-01 

2.71E-01 

2.04E-01 

29si 

8.99E-03 

6.03E-05 

8.61E-04 

6.05E-04 

3.87E-04 

2.49E-04 

30Si 

8.05E-03 

3.09E-05 

1.74E-03 

1.07E-03 

6.35E-04 

3.94E-04 

31p 

1.21E-03 

8.51E-05 

4.18E-04 

2.67E-04 

1.80E-04 

1.23E-04 

32s 

3.84E-02 

9.19E-02 

8.41E-02 

2.09E-01 

1.65E-01 

1.24E-01 

33s 

1.78E-04 

5.83E-05 

4.50E-04 

3.48E-04 

2.49E-04 

1.74E-04 
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Species 


Synthesized 

mass {Mq) 



Type II 
10~50 Mq 

W70 

W7 

Type la 
WDDl 

WDD2 

WDD3 


2.62E-03 

2.84E-06 

1.90E-03 

3.42E-03 

2.50E-03 

1.75E-03 

36s 

1.78E-06 

1.09E-11 

3.15E-07 

2.29E-07 

1.33E-07 

8.58E-08 

35 Cl 

l.OlE-04 

8.06E-06 

1.34E-04 

1.21E-04 

9.83E-05 

7.67E-05 

37C1 

1.88E-05 

5.36E-06 

3.98E-05 

4.23E-05 

3.36E-05 

2.52E-05 

36 A j- 

6.62E-03 

1.99E-02 

1.49E-02 

4.12E-02 

3.35E-02 

2.50E-02 

38 Ar 

1.37E-03 

5.93E-07 

1.06E-03 

1.71E-03 

1.45E-03 

1.12E-03 

40 A J- 

2.27E-08 

1.14E-12 

1.26E-08 

8.16E-09 

6.08E-09 

4.53E-09 

39k 

6.23E-05 

1.82E-06 

8.52E-05 

1.05E-04 

9.00E-05 

7.00E-05 

41k 

5.07E-06 

5.33E-07 

7.44E-06 

8.39E-06 

7.12E-06 

5.53E-06 

^OCa 

5.77E-03 

1.95E-02 

1.23E-02 

4.07E-02 

3.45E-02 

2.58E-02 

42Ca 

4.23E-05 

1.55E-08 

3.52E-05 

4.69E-05 

4.06E-05 

3.24E-05 

^^Ca 

1.08E-06 

6.81E-08 

1.03E-07 

6.97E-08 

6.31E-08 

4.73E-08 

^^Ca 

5.53E-05 

1.55E-05 

8.86E-06 

2.06E-05 

2.07E-05 

1.42E-05 

46 Ca 

1.43E-10 

5.88E-11 

1.99E-09 

6.17E-08 

7.18E-08 

3.00E-08 

48 Ca 

5.33E-14 

5.96E-12 

7.10E-12 

4.21E-06 

4.41E-06 

3.62E-06 

45Sc 

2.29E-07 

3.92E-08 

2.47E-07 

3.54E-07 

3.23E-07 

2.51E-07 

46 Ti 

7.48E-06 

7.87E-07 

1.71E-05 

2.06E-05 

1.76E-05 

1.44E-05 

47Ti 

2.11E-06 

7.49E-07 

6.04E-07 

1.21E-06 

1.24E-06 

7.40E-07 

48 Ti 

1.16E-04 

3.21E-04 

2.03E-04 

8.17E-04 

8.53E-04 

6.76E-04 

49 Ti 

5.98E-06 

1.64E-06 

1.69E-05 

6.55E-05 

6.71E-05 

5.24E-05 

50Ti 

3.81E-10 

1.14E-05 

1.26E-05 

4.57E-04 

4.60E-04 

4.80E-04 

50y 

7.25E-10 

5.66E-09 

8.28E-09 

5.49E-08 

5.69E-08 

4.52E-08 

Sly 

l.OOE-05 

2.10E-05 

5.15E-05 

2.92E-04 

3.18E-04 

2.10E-04 

50 Cr 

4.64E-05 

7.61E-05 

2.71E-04 

5.86E-04 

5.24E-04 

3.98E-04 

52 Cr 

1.15E-03 

6.63E-03 

5.15E-03 

1.75E-02 

2.01E-02 

1.67E-02 

53 Cr 

1.19E-04 

4.65E-04 

7.85E-04 

2.02E-03 

2.26E-03 

1.87E-03 

54 Cr 

2.33E-08 

1.79E-04 

1.90E-04 

2.03E-03 

2.03E-03 

2.08E-03 

55 Mn 

3.86E-04 

6.27E-03 

8.23E-03 

1.57E-02 

1.88E-02 

1.44E-02 

64Fe 

3.62E-03 

8.18E-02 

1.04E-01 

7.26E-02 

7.08E-02 

6.29E-02 

56 Fe 

8.44E-02 

6.72E-01 

6.13E-01 

4.56E-01 

6.15E-01 

7.62E-01 

s^Fe 

2.72E-03 

1.98E-02 

2.55E-02 

1.03E-02 

1.39E-02 

1.94E-02 

58 Fe 

7.22E-09 

9.34E-04 

9.63E-04 

4.09E-03 

4.06E-03 

4.10E-03 

69 Co 

7.27E-05 

6.28E-08 

1.02E-03 

6.83E-04 

8.60E-04 

1.12E-03 

58Ni 

3.63E-03 

9.67E-02 

1.28E-01 

3.01E-02 

3.34E-02 

4.62E-02 

eoNi 

1.75E-03 

1.43E-02 

1.05E-02 

4.03E-03 

4.15E-03 

4.78E-03 

61 Ni 

8.35E-05 

2.30E-04 

2.51E-04 

7.93E-05 

9.23E-05 

1.17E-04 

62 Ni 

5.09E-04 

1.37E-03 

2.66E-03 

1.26E-03 

1.36E-03 

1.88E-03 

64Ni 

3.20E-14 

1.22E-06 

1.31E-06 

3.45E-04 

3.35E-04 

3.61E-04 

63 Cu 

8.37E-07 

2.86E-06 

1.79E-06 

3.91E-05 

4.25E-05 

3.49E-05 

66CU 

4.07E-07 

9.60E-07 

6.83E-07 

1.48E-05 

1.59E-05 

l.OOE-05 

64 Zn 

1.03E-05 

1.05E-04 

1.22E-05 

3.54E-07 

6.97E-07 

2.56E-06 

66 Zn 

8.61E-06 

8.65E-06 

2.12E-05 

3.15E-05 

3.18E-05 

3.74E-05 

67Zn 

1.52E-08 

3.96E-09 

1.34E-08 

3.55E-06 

3.90E-06 

1.98E-06 

68 Zn 

3.92E-09 

3.29E-09 

1.02E-08 

4.38E-07 

4.94E-07 

2.81E-07 





